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The majority electron density as a function 
of the Fermi energy is calculated in zinc 
blende, n-type GaSb for donor densities 
between 10 cm" and 10 cm" . These 
calculations solve the charge neutrality 
equation self-consistently for a four-band 
model (three conduction sub-bands at T, L, 
and X and one equivalent valence band at 
r ) of GaSb. Our calculations assume par- 
abolic densities of states and thus do not 
treat the density-of-states modifications 
due to high concentrations of dopants, 
many body effects, and non-parabolicity of 
the bands. Even with these assumptions, 



the results are important for interpreting 
optical measurements such as Raman 
measurements that are proposed as a non- 
destructive method for wafer acceptance 
tests. 
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1. Introduction 

Most interpretations of optical measurements on 
compound semiconductors such as GaSb require phys- 
ical models and associated input parameters that 
describe how carrier densities vary with dopant concen- 
trations and measured Fermi energies. In this paper, we 
report on a method that gives closed form analytic 
expressions for the carrier densities in the conduction 
sub-bands for GaSb at room temperature. The method 
is based on an iterative and self-consistent solution of 
the charge neutrality equation with full Fermi-Dirac 
statistics for the carriers at finite temperature and on the 
use of statistical analyses to give analytic expressions 
that represent the calculated data sets. 

The method reported here is related to earlier work 
on n-type GaAs presented in reference [1]. Reference 
[1] gives the results predicted by an effective two-band 
model, one equivalent conduction band and one equiv- 
alent valence band at F, that includes the densities of 
states modifications due to high concentrations of 



dopants and due to many-body effects associated with 
carrier-carrier interactions. The method given below for 
GaSb is a four-band model. But, because of computa- 
tional limitations, it does not include the densities of 
states modifications due to high concentrations of 
dopants and due to many-body effects. 



2. Theory 

The electron n and hole concentrations h in units of 
cm"^ at thermal equilibrium are given, respectively, by 



« = J /„ (E)p^(E)dE and h=j[l- /,( E)] p„( £)d E, 



(1) 



where /o(£) = {l+exp[(£:-£:F)/^B7]}"' is the Fermi- 
Dirac distribution function, Ep is the Fermi energy in 
eV, pj(E) and pXE) are, respectively, the electron den- 
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sity of states for the conduction band and the hole den- 
sity of states for the valence band, ^b is the Boltzmann 
constant, and T is the temperature in kelvins. The cal- 
culations incorporate the Thomas-Fermi expression for 
the screening radius, 






VpSE)-p,{E)\^E, 



(2) 



and the charge neutrality condition 

N^ = n-h, 



(3) 



to compute self-consistently the Fermi energy E^ and 
the screening radius r^ for given values of the ionized 
dopant concentration A^i and temperature T, The static 
dielectric constant is £ and the permittivity of free space 
is Sq, The ionized dopant concentration is positive for n- 
type material (donor ions) and negative for p-type 
material (acceptor ions). The results reported here are 
for uncompensated n-type material. The results for the 
screening radius r^ are not reported here because they 
are not needed to extract carrier concentrations from 
Raman measurements. 

In this paper, we use the four-band model that has 
three conduction sub-bands centered at the F, L, and X 
symmetry points in the Brillouin zone and one equiva- 
lent valence band centered at the F symmetry point. We 
do not include the detailed nonparabolicity of the GaSb 
energy bands at F. Unlike GaAs, the GaSb conduction 
F, L, and X sub-band masses and energy spacings are 
such that for donor densities of technological interest, 
the conduction sub-band at L is the one that is most 
populated. The non-parabolicity of the conduction F 
sub-band in GaAs is discussed in Ref [2]. If we were 
to use the Kane three level k- p model [2], which does 
not include the conduction sub-bands at L and X, we 
would be able to include the non-parabolicity of the 
conduction F sub-band. However, because the conduc- 
tion F sub-band band in GaSb is not the dominant band 
for determining the Fermi energy its non-parabolicity 
correction may not have a significant effect on the 
results given below and may lie within the uncertainties 
associated with the band masses quoted in the literature 
for GaSb. 

The heavy hole mass m^ and light hole mass wix^ for 
the two degenerate sub-bands at the top of the valence 
band are combined to give an effective mass 



'^vr=('^hh +'^ih 7 



3/2\2/3 



(4) 



for the valence topmost sub-band. The values of these 
parameters are given in Table 1 . 

The zero of energy is at the bottom of the conduction 
F sub-band. The bottoms of the conduction L and X 
sub-bands are, respectively, at ^'^l and E^, The top of 
the degenerate valence F sub-band is at -Eq, where Eq 
is the bandgap of GaSb. The split-off valence sub-band 
at F due to spin-orbit coupling and the non-parabolici- 
ty factor of the conduction F sub-band are neglected. 
The probabilities for typical carriers in equilibrium to 
occupy appreciably these states are low. This means 
that the Fermi energies should be sufficiently above the 
valence sub-band maximum at F. Placing exact limits 
on the Fermi energies for which the four-band model is 
valid would be tenuous, because knowledge of how the 
various sub-bands move relative to one another due to 
the dopant concentrations considered here and due to 
many body effects is not adequate. 

The general expression [3] for the temperature 
dependence of conduction sub-band minima relative to 
the top of the valence band at F is 



E, = E,,-[A[r%T+B^-\ 



(5) 



in units of eV, where / = F, L, or X. The values for the 
coefficients E^^, A^, and B^ are listed in Table 2. 

The general expression for the parabolic densities of 
states for electrons and holes per band extrema and per 
spin direction is given by 



p(E) = 



NA7iV4e 



(8;r')(^'/2m*mo)''' 



(6) 



where A^^ is the number of equivalent ellipsoids in the 
first Brillouin zone, the volume of the unit cell is 
V= a^, aL is the lattice constant, m* is one of the effec- 
tive masses listed in Table 1 for the appropriate band 
extrema, and m^ is the free electron mass. Because eight 
permutations of the wave vector in the (HI) direction 
exist, there are eight L sub-band ellipsoids with centers 
located near the boundary of the first Brillouin zone. 
Also, because six permutations of the wave vector in 
the (100) direction exist, there are six X sub-band ellip- 
soids with centers located near the boundary of the first 
Brillouin zone. Since about half of each ellipsoid is in 
the neighboring zone, the number of equivalent sub- 
bands A^cL for the L sub-band is four and the number of 
equivalent sub-bands A^^x for the X sub-band is three. 

In terms of a four-band model for room temperature 
n-type GaSb, the total density of states pJ^E^ for the 
majority carrier electrons in n-type GaSb then becomes 
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Table 1. Input parameters for intrinsic zinc blende GaSb at 300 K. The energies of the extrema of the conduction and valence sub-bands are ref- 
erenced to the bottom of the conduction sub-band at the F symmetry point in the Brillouin zone of the reciprocal lattice space. The mass of the 
free electron is mq. These data are from Ref [3] 



Parameter 



Symbol 



Value 



Units 



Lattice constant 

Dielectric constant in vacuum 

Static dielectric constant 

Bandgap 

Bottom of the conduction L sub-band 

Bottom of the conduction X sub-band 

Top of the degenerate valence F sub-band 

Spin-orbit splitting 

Top of the split-off (spin-orbit splitting) valence F sub-band 

Effective mass of conduction F sub-band 

Transverse L sub-band mass 

Longitudinal L sub-band mass 

Effective mass of conduction L sub-band 

Transverse X sub-band mass 

Longitudinal X sub-band mass 

Effective mass of conduction X sub-band 

Light hole mass of degenerate valence F sub-band 

Heavy hole mass of degenerate valence F sub-band 

Effective mass of degenerate valence F sub-band 

Splitoff band mass of the valence sub-band at F 

Number of equivalent conduction L sub-bands 

Number of equivalent conduction X sub-bands 



«L 




6.09593x10"^ 


cm 


£ 




8.854x10"^^ 


F/m 


^0 




15.7 




^G^l^vrl 




0.726 


eV 


^cL 




0.084 


eV 


^vX 




0.31 


eV 


-^vr 




-0.726 


eV 


^so 




0.80 


eV 


-^sor"-^vr- 


-^so 


-1.526 


eV 


"^cr 




0.041 


niQ 


"^tL 




0.11 


niQ 


"^IL 




0.95 


niQ 


"^cL ^ (^IL "^tL 


2.1/3 


0.226 


niQ 


m^x 




0.22 


niQ 


mjx 




1.51 


niQ 


"^cX " ("^IX "^tX ) 


0.418 


rriQ 


m^ 




0.05 


niQ 


"^hh 




0.4 


niQ 


m^Y 




0.41 


niQ 


"^so 




0.14 


niQ 






4 
3 





p,(£) = p,r(£) + pj£) + p^(£), 



(7) 



where pci{E), paiE), and p^AE) are the sub-band den- 
sities of states for the conduction F, L, and X sub-bands 
with effective masses of mcr, rn^i., and m^^, respectively. 
The density of states for the minority carrier holes is 

PXE) = pAE) (8) 

with an effective mass of m^r- 

3. Results 

Tables 1 and 2 contain the input parameters for the 
calculations of the Fermi energy as a function of the 
dopant donor density. We solve self-consistently, by 
means of an iterative procedure, Eq. (3) with Eqs. (6), 
(7) and (8). The independent variable is the temperature 
T. The Fermi energy is varied for a given temperature 
until Eq. (3) is satisfied. Figure 1 presents the calculat- 
ed data graphically for 28 values of donor densities 
between 10^^ cm"^ and 10^^ cm"l Figure 2 gives the 
electron densities in the conduction sub-bands at F and 
L and the total electron density as functions of the 
Fermi energy. Figure 2 does not show the electron den- 
sity in the conduction sub-band at X, because it is less 



Table 2. CoeflRcients for the temperature dependence of the conduc- 
tion band extrema that are used in Eq. (5). These data are from Ref 
[3] 



Parameter 


Symbol 


Value 


Units 


r sub-band 


^ro 


0.813 


eV 


r sub-band 


Ay 


3.78 X 10^ 


eV/K 


r sub-band 


By 


94. 


K 


L sub-band 


^LO 


0.902 


eV 


L sub-band 


Al 


3.97x10"^ 


eV/K 


L sub-band 


Bl 


94. 


K 


X sub-band 


^xo 


1.142 


eV 


X sub-band 


Ax 


4.75 X 10^ 


eV/K 


X sub-band 


Bx 


94. 


K 



than 10"^ times the total electron density Because 
m^Y ^^ '^cL and E^Y is much closer to E^y than it is to E^x> 
the electron density in the conduction L sub-band 
exceeds the electron density in the conduction F sub- 
band at room temperature. The solid curve in Fig. 2 is 
the same curve as given in Fig. 1. Figure 2 shows that 
the majority of electrons is in the conduction L sub- 
band and that the density of electrons in the L sub-band 
approaches the total density of electrons as the donor 
density approaches 10^^ cm"-^. Hence, even though 
GaSb is intrinsically a direct semiconductor, the results 
from Fig. 2 suggest that electrons for n-type GaSb in 
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Fig. 1. The calculated Fermi energy for n-type GaSb at 300 K as a 
function of the donor density. The Fermi energy is relative to the 
majority conduction band edge at the F symmetry point in the first 
Brillouin zone. 
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Fig. 2. The calculated electron densities in the conduction sub-bands 
at r and L and the total electron density as functions of the Fermi 
energy. Tbe Fermi energy is relative to the majority conduction band 
edge at the r symmetry point in the first Brillouin zone. 



the vicinity of the Fermi surface will behave as though 
they have many characteristics of electrons in an indi- 
rect semiconductor 

For illustrative purposes, we give here only the 
results for fitting the logarithm to the base 1 of the 
total electron density n in units of cm~^ to a polynomial 
in £'f, namely 



log JO (« cm-'): 



■■ a^^ + a^Ep + a^Ep +. 



. + ^^.... (9) 



The analytic fits for the electron densities in the F, L, 
and X sub-bands are available by sending an email to 
herbert.bennett@nist.gov. During the fitting analyses, 
we rely substantially on graphics and keep the number 
of fitting parameters to a minimum, subject to the con- 
straint that the residual standard deviation S.^^ is accept- 
ably small, i.e., 5'res<0.01. The standard deviation is a 
measure of the "average" error in a fitted model and 
thereby is a metric for assessing the quality of the fit. A 
smaller S^ indicates a better fit. The residual standard 
deviation for a model Y^=f{Z) is 



E(V 



-Y/f/(N- 



P)l 



(10) 



7=1 



where Yj are the calculated data values, the Yf are the 
predicted values from the fitted model, A^ is the total 
number of data points (here A^= 28), and P is the total 
number of parameters to be fitted in the model. We use 
the NIST-developed DATAPLOT [4] software for both 
the exploratory graphics and for the statistical analyses. 
Table 3 gives the four fitting parameters for the cubic 
/ = 3 polynomial fit to the logio(^) ^s shown in Eq. (9) 
and the associated residual standard deviation S,^^ = 
0.0066. In general, the values ofS,^^ decrease monoton- 
ically with increasing number / of terms in the polyno- 
mial given in Eq. (9). But, care must be taken to avoid 
fitting noise in data sets. The general guideline for 



Table 3. The four fitting parameters for a cubic polynomial fit of the theoretical calculations for the total 
electron density in n-type, zinc blende GaSb at 300 K as a function of the Fermi energy relative to the bot- 
tom of the conduction T sub-band. The ratio is the estimated value divided by the estimated standard devia- 
tion. The residual standard deviation is S^^ = 0.0066 



Fitting parameter 



Estimated value 



Estimated standard 
deviation 



Units 



Ratio 



aQ 


17.7504 


0.1774x10" 


a' 


15.6775 


0.5416x10" 


a^ 


-11.4745 


0.4723 X 10" 


a' 


^1.3848 


0.8535 



eV- 



1.001 X 10' 
2.895 X 10^ 
-2.43 X 10^ 
-4.849 X lO' 
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many data sets is that when the absolute value of the 
ratio R of the estimated parameter value divided by its 
estimated standard deviation is less than about 2, then 
the rate of decrease in S.^^ with increasing / tends to 
decrease. For the data given in Fig. 1, when / = 2 or 
P = 3, 5'res = 0.0644; when 1 = 3 or P = 4, S,,, = 0.0066; 
and when 1 = 4 or P = 5, 5*^^3 = 0.0063. Because the 
change in values of S.^^ between 1 = 3 and / = 4 is not 
significant, we use the fitting parameters for the cubic 
1 = 3 case in this paper. Also, the ratio R for the param- 
eter a4 when 1 = 4 is -1.917, and such a value for R 
means that proceeding with higher / values probably is 
not warranted. A figure that compares the calculated 
total electron density as a function of the Fermi energy 
with the fitting results from Eq. (9) for a cubic polyno- 
mial is not given because the two curves lie on top of 
one another within the line widths of each curve. Fits to 
the calculated electron densities for each of the conduc- 
tion sub-bands are available from the author upon 
request. Also, since the screening radii for the carriers 
from Eq. (2) are not needed when interpreting the pro- 
posed measurements considered here, the correspon- 
ding screening radii are not presented in this paper. 



4. Conclusions 

The results from Sec. 3 are consistent with the find- 
ings of experimental work reported in the literature 
such as Refs. [5] and [6]. Interpreting experiments for 
GaSb requires at least a three-band model and under 
some conditions may require a four-band model. And 
finally, even though GaSb is intrinsically a direct semi- 
conductor, our results show that electrons for n-type 
GaSb in the vicinity of the Fermi surface will have 
some characteristics that are similar to those for elec- 
trons in an indirect semiconductor. 
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